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We consider Clifford Quantum Cellular Automata (CQCAs) and their time evolution. CQCAs 
are an especially simple type of Quantum Cellular Automata, yet they show complex asymptotics 
and can even be a basic ingredient for universal quantum computation. In this work we study 
the time evolution of different classes of CQCAs. We distinguish between periodic CQCAs, fractal 
CQCAs and CQCAs with gliders. We then identify invariant states and study convergence prop- 
erties of classes of states, like quasifree and stabilizer states. Finally we consider the generation of 
entanglement analytically and numerically for stabilizer and quasifree states. 
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I. INTRODUCTION 



Quantum cellular automata (QCAs), i.e., reversible quantum systems which are discrete both in 
time and in space and exhibit strictly finite propagation, have recently come under study from 
different directions. On the one hand, they serve as one of the computational paradigms for quantum 
computation, and it has been shown that a certain one-dimensional QCA with twelve states per 
cell can efficiently simulate all quantum computers @j. On the practical side, QCAs are a direct 
axiomatization of the kind of quantum simulator in optical lattices, which are under construction in 
many labs at the moment [H H[. Related to this, they can be seen as a paradigm of quantum lattice 
systems, in which the consequences of locality, assumed in the idealized pure form of strictly finite 
propagation, can be explored directly. Due to the famous Lieb-Robinson bounds 043 this feature is 
also present in continuous time models albeit in an approximate form. 

In all these settings, the time asymptotics for the iteration of the QCA is of interest, and displays 
a curious dichotomy between a global and a local point of view. On the one hand we are assuming 
reversibility, so the global evolution is an automorphism taking pure states to pure states. If we split 
the system into two subsystems, e.g., a right half chain and a left half chain, then we expect the 
QCA to generate entanglement from any initial product state. There is a simple upper bound (see 
Section TlV A[) showing that the entanglement growth is at most linear, and we find indeed that for 
the automata studied in this paper this is the typical behavior. However, utilizing this entanglement 
requires the control of larger and larger regions. So from a local perspective, i.e., when we only 
consider the restriction of the state to a finite region, we will not see this increase. In fact, a typical 
behavior of the local restrictions is the convergence to the maximally mixed state. That for large times 
the state seems globally pure and locally completely mixed is no contradiction: it merely reflects the 
fact that any local system becomes maximally entangled with its environment. 

Stationary states are in some sense the final result of an asymptotic evolution. Again the local 
analysis makes it clear that the totally mixed state is invariant for any reversible cellular automaton. 
In general there may be many more invariant states, among them some, which are not only ergodic 
(i.e., extremal in the set of all translation invariant states), but even pure (i.e., extremal in the set of 
all states). For a special class we consider here, we exhibit a rich set of such states. 

The very fact that QCAs can serve as a universal computational model suggests that asymptotic 
questions cannot be easily answered in full generality. Therefore, in this paper we look at a subclass 
of QCAs, the Clifford Quantum Cellular Automata (CQCAs). In this case much of the essential 
information can be obtained by studying a classical cellular automaton, which even turns out to be of 
a linear type. Consequently, we can answer some questions exhaustively. The drawback is, of course, 
that universal computation is not possible in this class. However, we believe that some typical features 
of QCA asymptotics can be studied in this theoretical laboratory . 

Our paper is organized as follows. We begin with a short introduction to CQCAs and their repre- 
sentation as 2 x 2-matrices in Section III A[ Then we derive a general theory of the time evolution of 
CQCAs in Section III Bl We show that CQCAs can be divided in three major classes by their time 
evolution: periodic automata (Section III B 2[) . automata which act as lattice translations on special 
observables we call gliders (Section III B and fractal automata (Section III B 3[) . whose space-time 
picture shows self-similarity on large scales. The class of a CQCA is determined by the trace of its 
matrix. A constant trace means the automaton is periodic, a trace of the form u~ n + u n indicates 
gliders that move n steps on the lattice each time step. All other CQCAs show fractal behavior. We 
prove that all automata with gliders and n = I are equivalent in the sense that they can be trans- 
formed into each other by conjugation with other CQCAs. This turns out not to be true for glider 
CQCAs with n> 1. 
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Using the results for the observable asymptotics, in Section IIIII convergence and invariance of 
translation-invariant states are analyzed. For periodic CQCAs the construction of mixed invariant 
states is straightforward (Section IIII Aj) . There are also pure stabilizer states that are invariant with 
respect to periodic CQCAs (Section IIII C ip . Non-periodic CQCAs can not leave stabilizer states 
invariant. In fact, for fractal CQCAs the only known invariant state is the tracial state. For glider 
CQCAs invariant states exist, and are constructed as the limit of the time evolution of initial product 
states in Section HUH Finally we consider n = 1 glider automata and states which are quasifree on 
the canonical anticommutation relations (CAR) algebra (Section llll D[) . We employ the Araki-Jordan- 
Wigner transformation (Section [ill D II) to transfer the glider CQCA to the CAR algebra and study 
its action as a Bogoliubov transformation on quasifree states; we study both convergence (Section 
IIII D 3|) and invariance (Section lniD2j) . 

In the last Section (Section IIV[) the entanglement generation of CQCAs is considered. First we 
derive a general linear upper bound for entanglement generation of QCAs. We then prove that in the 
translation-invariant case the bound is more restrictive and can be saturated by CQCAs acting on 
initially pure stabilizer states. For CQCAs acting on stabilizer states the asymptotic entanglement 
generation rate is governed by the order of the trace polynomial (Section IIVB[) . We then examine 
a set of quasifree states interpolating between a product state and a state which is invariant with 
respect to the standard glider CQCA (Section IIVC|) . It is shown numerically, that the entanglement 
generation is linear also in this case, but with a slope that can be arbitrarily small. 

II. TIME EVOLUTION OF THE OBSERVABLE ALGEBRA 

CQCAs show a variety of different time evolutions. In this section we develop criteria to predict 
the time evolution by characteristics of the matrix of the corresponding symplectic cellular automaton 
(SCA). 

A. A short introduction to CQCAs 

Clifford Quantum Cellular Automata are a special class of Quantum Cellular Automata first de- 
scribed in [9(. As the name indicates, they are QCAs that use the Clifford group operations. They 
can be defined for arbitrary lattice dimensions and prime cell dimensions. Here we will only consider 
the case of a one dimensional infinite lattice and cell dimension two (qubits). Thus we deal with 
an ordinary spin-chain. QCAs are translation-invariant operations which we define on a quasi-local 
observable algebra Reversible QCAs are automorphisms. So in our case CQCAs are translation- 
invariant automorphisms of the spin-i chain observable algebra 21 = <8>^_ 00 2ti, where 21^ = M 2 . In 
explicit this means that a CQCA T commutes with the lattice translation r and leaves the product 
structure of observables invariant: 

T(AB) = T(A)T(B), VA, B e 21 

As shown in [l| a QCA is fully specified by its local transition rule To, i.e., the picture of the one-site 
observable algebra. In our case the Pauli matrices form a basis of this algebra, so their pictures specify 
the whole CQCA. There is one important restriction on the set of possible images. They have to fulfill 
the commutation relations of the original Pauli matrices on the same site and also on all other sites. 
As a QCA in general has a propagation on the lattice, one cell observables are mapped to observables 
on a neighborhood Af of the site. Thus neighboring one-site observables may overlap after one time 
step but still have to commute. This imposes conditions on the local rule. CQCAs are defined as 
follows: 

Definition II. 1. A Clifford Quantum Cellular Automaton T is an automorphism of the quasi-local 
observable algebra of the infinite spin chain that maps tensor products of Pauli matrices to multiples 
of tensor products of Pauli matrices and commutes with the lattice translation r. 

We now want to find a classical description of the CQCA. It is well known that Clifford operations 
can be simulated efficiently by a classical computer. Therefore it is not surprising that an efficient 
classical description of CQCAs exists. This description was introduced in [9(. We will only give a 
short overview of the topic, for proofs and details we refer to the literature. 
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We use the (finite) tensor products of Pauli matrices as a basis of the observable algebra. Every 
CQCA T maps those tensor products to multiples of tensor products of Pauli matrices. The factor 
can only be a complex phase, that can be fixed uniquely by the phase for single cell observables (single 
Pauli matrices). We can thus describe the action of the CQCA T on Pauli matrices by a classical 
cellular automaton t acting on their labels "1, 2, 3". We could keep track of the phase separately, but 
for our analysis this is unnecessary. 

Mathematically we describe this correspondence as follows: The Pauli matrices correspond to a 
Weyl system over a discrete phase space where all operations are carried out modulus two. For one 
site we have 

CT1 = X = W(1,0), C72 = y = iw(l,l), C73 = Z = w(0, 1), (70 =1 =w(0,0). 
Tensor products of these operators are constructed via 

w(0 = ®w(e(a:)) 

where £ = (£+ , £_ ) is a tuple of two binary strings which differ from on only finitely many places and 
is its value at position x, e.g. (J). Thus the tensor product is well defined. Before we continue 
with the mathematical definition, we want to illustrate the classical description by a simple example: 

Example II. 2. We define our CQCA on the observable algebra 21 of a spin chain by the rule 

Ti : 21, -> 2W+i, 
Tcrl = 4, 
Ta\ = al' 1 ® a\ ® a l 3 +1 . 

The image of oi follows from the product of the images of <j\ and 03 : 

Tcr| = -o^ 1 ®a\® al +1 . 

T has to be an automorphism to be a CQCA. To verify this we check if the commutation relations 
are preserved: 

[Ta[,Ta{] = [4,4] =0, 

[t4,t4) = wr 1 ® 4 ® <r^- l ,4~ 1 ® 4 ® 4 +1 ] = °> 

and 

[T4 , Ta{] = ®*\®4 +1 ,4} = i+ j, 
{T4, Ta{} = ®a\® 4 +1 ,4} = i = j . 

This automaton will be used extensively in the following parts of the paper, so we give it the name 
G s . If we think of the CQCA G s as a classical automaton acting on the labels of the Pauli matrices 
we can illustrate the evolution (for one time step) of the observable a^ 1 ® 0% ® u\ as follows (the 
underlined labels are situated at the origin): 

3 1 3 

0.(3 21) - 9s (300) • ffs (020) • 9s (001) = ® 3 — 3 = (321). 

= 3 -«2 il 

We observe, that the observable only moves on the lattice under the action of the CQCA T. We call 
observables with this property gliders. Their existence can be observed easily, when we consider the 
space-time images of one-cell observables. o\ and 173 generate "checkerboards" of <j\ and (73 matrices. 
As o\ is mapped to 173 in the first step, the <7i-checkerboard is the same as the (73-checkerboard shifted 
one step in time. If we additionally shift in the space direction by one cell, the checkerboards are 
exactly the same up to two diagonals and thus cancel out as shown in Figure [TJ We thus produced a 
very simple observable on which the automaton acts as a translation, a basic glider. 

Another interesting property of this automaton is the fact that it maps the "all spins up" product 
state to a one-dimensional cluster-state, which is a one-dimensional version of the two-dimensional 
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FIG. 1: Glider of the example CQCA JT]). The blue boxes represent o\, the yellow ones 03. 



resource-state for the "One Way Quantum Computing" scheme by Raussendorf and Briegel 10]. It 
is also the basic ingredient of a scheme of "quantum computation via translation-invariant operations 
on a chain of qubits" by Raussendorf [111 ]. In a similar way, the update rule G s (but with <j\ and 
03 exchanged) has appeared as time-evolution of spin chains implemented by a Hamiltonian that is 
subjected to periodic quenches [lU, [l3j], and it has even been realized experimentally in an NMR- 
System pi]. 

In the phase space picture we can describe the automaton by a 2 x 2-matrix with polynomial entries. 
In phase space our CQCA-rule reads 

'i\ f o' 

9 S 



oj y 1 J 

0-(;i!)- 

D -(;;:)■ 

Now we transform the binary strings to Laurent polynomials by indicating the position by a multipli- 
cation with a variable u and add all terms from the different positions to a Laurent polynomial: 1 

3 ■ G 

0\ I 1 

1/ \u- 1 +u 



') = ( 1 



a 
a 
a 

We arrange the images of (J) and L) in a 2 x 2-matrix 

(1) 

The image under G s of an arbitrary tensor product of Pauli matrices is now determined by the 
multiplication of the corresponding vector of polynomials by the matrix representation g s of G s . We 
will later argue that this works for all CQCA. 

Now we come back to the mathematical definition of CQCAs: The Weyl operators fulfill the relation 

w« + »() = (-l)^wK)w(,) 




1 More formally, we perform an algebraic Fourier transformation 9] mapping the binary string £(x) to the Laurent 
polynomial = ExSZ^V- 
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and therefore the commutation relation 

w(Ow(77) = (-1)«+ i '--«-"+w(t7)w(£). 

In both cases terms of the type £+77- are scalar products where the addition is carried out modulo 2. 
The arguments of the Weyl operators are elements of a vector space over the finite field Z2 which we 
call phase space and thus commute. Of course the corresponding Weyl operators do not necessarily 
commute, but they always commute or anticommute. Their commutation relations are encoded in 
the symplectic form <r(£, if) = n) — /3(ry, £) = £+77- — £-?7+ £ Z2. As an automorphism the CQCA 
leaves the commutation relations invariant. A representation of the CQCA on phase space therefore 
has to leave the symplectic form invariant. Such a translation- invariant symplectic map is called 
symplectic cellular automaton (SCA). We can find a SCA and an appropriate phase function A(£) for 
every CQCA. 

Proposition II. 3 ([9]). Let T be a CQCA. Then we can write 

r[w(0] = A(Ow(t0 (2) 

with a symplectic cellular automaton t and a translation invariant phase function A(£) which fulfills 

A(£ + V ) = A(C)A(77)(-l)^^^ (t «^ 

as well as |A(£)| = 1 V£. Furthermore, A(£) is uniquely determined for all £ by the choice of X on one 
site. 

In the following analysis of CQCAs we neglect the global phase and consider the symplectic cellular 
automata only. As we can always find appropriate phase functions all results for SCAs translate to the 
world of CQCAs directly. We have already seen in Example 111.21 that there exists a very convenient 
representation of CQCAs as 2 x 2-matrices with polynomial entries. 

Definition II. 4. V is the ring of Laurent polynomials over Z2. TZ is the subring ofV, which consists 
of all polynomials, which are reflection invariant with center u = 0. 

Theorem II. 5. Every CQCA T is represented up to a phase by a unique 2 x 2-matrix t with entries 
from V . Such a matrix represents a CQCA if and only if 

• det(t) = u 2a , a e Z; 

• all entries are symmetric polynomials centered around the same (but arbitrary) lattice point a; 

• the entries of both column vectors, which are the pictures of (1, 0) and (0, 1), are coprime. 

Proof. We will only give a sketch of the proof. For details see [9J. The connection between CQCAs 
and SCAs was already established in Proposition 111.31 What remains to show is that SCAs are linear 
transformations over V 2 . The application of a SCA to a vector in phase space can be described as 
the multiplication of this vector with a matrix representing the SCA from the left. The product is 
the convolution of binary strings. By the algebraic Fourier transform f(u) — X^ez f( x ) uX i which 
maps the vectors of binary strings onto vectors with entries from the ring V of Laurent polynomials 
over the finite field Z2, the convolution of strings translates into the multiplication of polynomials. 
Thus in this picture the application of the SCA to a phase space vector is just a common matrix 
multiplication, where all operations are carried out modulo 2. If we translate the symplectic form and 
the condition that it has to be left invariant to the polynomial picture we retrieve the above conditions 
on the matrix t. □ 

We can further simplify these statements by only considering automata centered around 0. The 
lattice translation r is a SCA which by definition commutes with all other SCAs. Its determinant 
is given by det(r) = u 2 . Therefore every SCA can be written as the product of a lattice translation 
and an automaton centered around which has determinant one. We call these automata centered 
symplectic cellular automata (CSCAs) and in the following sections we point our focus to those. 

CSCAs and CQCAs each form a group. This group is generated by a countably infinite set of basis 
automata. The CSC A form the group T = SL(2,72.), which is the group of all 2 x 2-matrices with 
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determinant 1 over the ring of centered Laurent polynomials with binary coefficients. The group Tq 
is the group of local automata. Their generators are 




f 



Additionally we have the shear transformations 



G n :=\ 1 ° ] , n G N. 
\u n + u- n 1 J 

which complete the set of generators of T. For proofs see 

B. Classification of Clifford quantum cellular automata 

The time evolution of a CQCA is determined up to a phase by the powers of the matrix of the 
corresponding CSCA. We will only consider the evolution of single cell observables, as any other 
observable can be represented by products and sums of these. This product structure is invariant 
under the action of the automaton, because it is an automorphism of the observable algebra. This 
means that when we discuss the time evolution of a CQCA T, we will consider the action of powers 
t n of the matrix of the CSCA on phase space vectors which only contain constants. For example the 
image of o\ — w(l, 0) after n time steps of T is given by the first column vector of t n (and a global 
phase). The matrix t does not always have an eigenvalue, because it is a matrix over a ring without 
multiplicative inverses for all elements. But for some of the automata eigenvalues do exist. These 
automata are called glider-automata, because on a special set of observables, the gliders, they act as 
lattice translations. 

We will prove that if the trace is a polynomial consisting of only two summands, i.e., it is of the 
form tr t — u~ n + u n , two eigenvectors exist and the automaton has gliders. If this does not hold, we 
can distinguish two cases. The trace can be either a constant, or an arbitrary symmetric polynomial. 
In the first case the automaton is periodic, in the second case it generates a time evolution which has 
fractal properties. As the case of periodic automata is not very interesting and the case of fractal 
automata will be covered in another paper fl5j . we focus on automata with gliders. We prove that all 
automata with gliders which move one step on the lattice at every timestep are equivalent. We also 
give an example to show that this is not the case for gliders that move more than one step. 

1. Automata with gliders 

We will first define our notion of a glider. Here we consider the case of qubits only, but with minor 
alterations all results of this section also hold for qudits with prime dimension. This extension and 
also all proofs which are omitted here can be found in [l|| [ljj ■ 

Definition II. 6. A glider is an observable, on which the CQCA acts as a lattice translation. In the 
Laurent polynomial picture a translation is a multiplication by u k , k G Z. 

We have already seen this behavior in Example 111.21 Now we determine the conditions a CQCA 
has to fulfill to have gliders. In general we can not diagonalize the matrices of the corresponding 
CSCA, because all our calculations are over a finite field and the entries are only polynomials in u 
which are centered palindromes. The polynomials are usually not invertible, so the equations which 
occur in the diagonalization cannot be solved mechanically. Furthermore, the diagonal matrix would 
not correspond to a CSCA as u n and u~ n are not centered palindromes. Hence we take a different 
approach. First, let us introduce some terms: We call a glider £ = (£+,£_) a minimal glider iff its 
two entries in phase space £+, £_ have no common non- invertible divisor 2 . The wedge-product of two 



The phase space vector of a minimal glider is maximal with respect to the notation introduced in Q. 
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phase space vectors shall be defined as £ A r\ = £+?7_ — ?7+£-- As we deal with qubits here, addition 
and multiplication of polynomials are carried out modulo two and £ A r\ = £+77- + ?7+£-- We define 
the involution of a polynomial p as the substitution of u by u^ 1 and denote it by p. Finally we will 
also need the following proposition: 



Proposition II. 7 ( [16l Il7j|). In the ring V of Laurent polynomials over the finite field Z 2 , the only 
invertible elements are monomials. 

Now we have all we need to state following theorems: 

Proposition II. 8. Given a CSCA t and a non-zero phase space vector £ with t£ = u n £, neH the 
following is true: 

1. £ fulfills t£ = it~ n £, thus it is a glider with the same speed but different direction as £. 

2. t is uniquely given by 



L 'U = — : (3) 

(4) 
(5) 

t 22 = «-^-+»"^- . (6) 



£A£ 

I- + 
£A£ 



5. tr(t) = vr n + tt n 

^. All gliders are multiples of 



or 





/ gcd(«' 


'+tn,t 12 ) \ 


1 


•"+tn / 


gcd(u T 


'+tu,ti 2 ) / 




tl2 \ 


/ gcd(«- 


"+tll,ti 2 ) \ 


V u ~ 


"+tll J 


gcd(u _ 


"+tll,ti 2 ) / 



(7) 



(8) 



Proof. 1. We use i£ = u™£ and take the involution on both sides, t consists of palindromes, thus 
i = t and we get t£ = 

2. We write = it™£ and i£ = u~ n £ component wise yielding the four equations 

(I) tuCf + tiaC- = 

(II) i 2i e+ + <22^- = U»£_, 

(I) tnCf + tue- = 

(II) t 2 l£f + t 2 2£- = U-"|L. 

Combining them in the right way we get 

hi (£ A £) = u B £+£_ + (9) 

«i2(fAC) = («" + «-")£+£+, (10) 

*2i(eA?) - + (11) 

* 22 (€ A |) = + (12) 

By assumption i is a CSCA, so the division by £ A £ gives a polynomial result and we have 
Equations §S§ to ^ . In Proposition III. 101 we will show which conditions £ has to fulfill for the 
division to be valid and therefore £ to be a glider. 
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3. 



tit = U 



(«A€) 



4. We now use (I) and (II) together with det t = 1 and tr t = u ™ + u" to derive the form of £. 
We get the equation 

e+(« n -tll) = Ml2. 

This equation for and , has still one free parameter. One particular solution for the 
equation is ^ art ' = f 12) ^ < J ,art '> — u n _ ^ obtain the minimal glider we have to divide 
these components of the particular solution by their greatest common divisor, and thus obtain 
(0. For © we do the same with (I) and (II). An arbitrary glider can be written as the glider 
defined by either or (|8]) times a Laurent polynomial in u. 

□ 

Remark II. 9. We could extend our definition of gliders to observables with polynomial eigenvalues 
A. These observables would be mapped to products of translates of themselves. We can show that 
this extension would not yield any new gliders: A CSCA t has to fulfill det t = 1. With (/) to (II) 
we get 

deU-(£A£) 2 = (tiit 22 + ti2*2i) ■ (£A£") 2 

= (A£+ f- + A£+£_) • (A£ + £_ + A£+£_) 

-((A + A) ((A + A) 

= AA + 2£+ + 

= AA(^A^) 2 

det t = AA = 1 

<S> A = u". 

The only possible solutions are A = u ±n , n £ N, because by Proposition III. 71 in only monomials 
have inverse elements. We know that £ is always a glider to — n so we only look at positive n. u° is 
excluded, because there is no propagation then. 

Proposition 11.10. A minimal £ £ V 2 is a glider for a CSCA t with eigenvalue u n , n £ Z\0 if and 
only if is a divisor of u~ n + u n . 

Proof. First let us assume, that £ is minimal, and £ A £ is divisor of u~ n + u n . For £ to be a glider 
with eigenvalue it™, — A£ has to hold. Therefore the division in Equations © to (JT2J) has to be 
valid. For the Equations (p~0|) and (|TT|) this is obviously true. For the other two equations we use a 
simple trick: 

tn (e a = + 

V v ' 

=0 

= + tr»(£ A £) 

Now it is apparent, that £ A £ also divides the right hand side of ([9]) if it is a divisor of u~ n + u n . For 
(fT2l) an analogous argument holds. 

Now let us show the converse. The Laurent polynomials form an Euclidean ring, which implies that 
the extended Euclidean algorithm is applicable [18j . Since £ is minimal, the greatest common divisor 
of £ + and £_ is 1, according to the extended Euclidean algorithm we can chose 77+ and 77 _ such that 

£+77- +£_tH- = gcd(e+, £-) = !■ 
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Then we have: 

( u - n + u n ) = ( u - n + u n ) ■ (£+77- + n-v+) ■ {i+fi- + i-n+) 

* . ' "> » ' 

= 1, as £ min. =1, as f min. 

= ( u - n + u n ) ■ (z+z+v-fi- + Z-Z+v+v- + Z+S-v-v+ + S-S-v+v+) 

= {u- n + u n ) C+Z+V-V- + (u- n + u n ) ■ £-t-ri + fj + 

S v ' 

=0 

S v ' 

=0 

= (u- n + u n ) £+£+ 77_77 + (u~ n + u n ) 77+77+ 

=*12(SAf) =*2l(fA£) 

+ («- n f+C- + « n £-£+ ) 77+77- + (u n £+ £_ + u-"C_e + ) 7?_r7+ 

V * ' " « ' 

=t2 2 (£Af) =in(?AC) 
+ U -"r7 + f?- (S A I) + ^"77-77+ (£ A £) 

which implies that £ A £ divides (« -n + u n ). □ 

We have shown in Proposition lII.8l that tr t = u~ n + u n is a necessary condition for t to have gliders. 
The following proposition shows that this condition is also sufficient. 

Proposition 11.11. A CSCA possesses gliders with eigenvalues A+ = u n and A_ = A+ = u~ n if and 
only if tr t = u~ n + u n . 

Proof. The "only if" part was already shown in Proposition III. 81 

We now assume that tr t = u~ n + u n and use this to evaluate the characteristic polynomial of t. 
Using det t — 1 we get 

A 2 + A • (iT" + u n ) + 1 = 0, 

which is solved by A± = u ±n , n e N. Thus the CSCA possesses gliders. □ 

Now that we know the conditions for the existence of gliders, we want to know how and when they 
can be connected. Consider an arbitrary CSCA t with gliders and a second CSCA b. If we transform 
t by conjugating with b we get t = b l tb which has by 

trt = tr(6 _: 4fo) = tr(lt) = tr t 

the same trace as t and thus is a glider automorphism, too. What is maybe more surprising is that 
the converse is also true for gliders of propagation speed one: any CSCA with one-step gliders is 
equivalent to the standard-glider CSCA g s ([T]) by the equivalence relation t = b~ 1 g s b. 

Theorem 11.12. Let £ = (£+, £_) 6 V 2 be minimal. Then the following three statements are 
equivalent: 

1. There is a CSCA t with t£ = u£. 

2. There is a CSCA b with b£ = Q . 

3. £A£ = (u- 1 +u). 

Proof. 3 1 has already been shown in Proposition III. 101 

1 => 2: We assume that 1 (and therefore also 3) is true and analyze the conditions this imposes on 
b: We know that b = b. We start with the assumption b£ = (*) and obtain 



bn(£A£) = £-+£-, (13) 

MCA^) - £++£+, (W) 

&2i(£A£) = (tT 1 ^ + u£_) , (15) 

622 (£A£) = + u£+). (16) 
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First we need to show that the matrix b actually exists, i.e., that all the right sides of the equations 
(fT3| to (fT6| can be divided by £ A £. This follows from the same argument as in Proposition III. 101 if 
tr t = u^ 1 + u, which is given by 1 3. For the matrix to be a CSCA the determinant has to be one. 
This is also true and can be shown by direct computation: 

det b = —^-2 ■ ((£_ + £_) ■ (u-^ + + < + ) + (u-^- + <-) • (£+ + £+)) 



i 



£A£ 



(u-'+u) ■ (£A£~) 



= 1 



This step only works for £ A £ = (it -1 + it) which corresponds to one step gliders. Later on we 
will consider gliders with higher propagation speed and give counterexamples to similar notions of 
equivalence for their automata. 
2 =*> 3: 

£ A £ = detft, ■(£ A £) = ft£ A &£ = Q A - u" 1 + u 

=i 

This completes our proof. □ 

We will now show, that for automata with higher propagation speed the gliders are not equivalent 
in the above sense. It is apparent from proposition III. 101 that for a fixed n > 1 there are gliders 
with different wedge products. These can not be transformed into each other, because the wedge 
product is invariant under transformation with CSCAs (see Theorem III. 121 part 2). Another way to 
see that there are different types of n-step glider automata is the fact, that we always have automata 
which are powers of one-step automata and also automata, whose roots are not CSCAs. These can 
not be transformed into each other. But even automata for gliders with the same wedge product 
can not always be connected by a third CSCA. To show this, it is sufficient to find two phase space 
vectors £ and r\ with £ A £ = rj A fj and £ A £ dividing (u~ n + u n ) for some n which can not be 
transformed into each other by a CSCA. We choose £ = ( u } u a) and rj = ( ^ u ). Their wedge product 
£ A £ = r] Afj — u~ 2 + u^ 1 +u + u 2 divides (u~ 3 + u 3 ). It is a valid wedge product for a 3-step glider. 
If an automaton b with br/ = £ existed, it would have to fulfill the equations 

&n?7+ + bx2V- = £+ = 1 
biif} + + &i2?7- = £+ = 1. 



From these we get 



M£A£) = V-+V- 
bn(u~ 2 + u~ x + u + u 2 ) = u~ 2 +u 2 

which can not be solved by any &n € 72.. 

2. Periodic automata 

CSCAs whose matrices have a trace independent of u show periodic behavior. 

Proposition 11.13. A CSCA t is periodic with period c + 2 if tit = c for c e Za. 

Proof. By the Cayley-Hamilton theorem we get t 2 = t ■ tr t + 1 = ct + 1 and thus t 2 = 1 for c = 
and t 3 = t 2 + t = 2t + 1 = 1 for c = 1. □ 

Proposition 11.14. Let t be a CSCA and £ G V 2 a non-zero phase space vector. If t£ = £ holds, 
then t is periodic with period two. 
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Proof. We use t 2 = t • tr t + 1 and i£ = £: 



t 2 i = Z 

{t-tit)z + tZ = z 

£(tri) = 

o tri = 

Thus t is of period two by Proposition III. 13l □ 



3. Fractal automata 



All CSCAs that are neither periodic nor have gliders show a fractal time evolution. Fractal means, 
that the graph of the spacetime evolution of one cell observables is self similar in the limit of infinitely 
many timesteps. We will cover this type of CSCAs in detail in a future publication Here we only 
give an example to illustrate the self similarity. The evolution of the automaton 



/ = 



M" 1 + 1 +U 1 

1 



(17) 



is shown in Figure [2j 




FIG. 2: Time evolution of a fractal CSCA (|17[) . Time increases upwards. The different colors illustrate the 
different Pauli matrices. 



Nevertheless, we will state one short lemma, that we will need later on to prove the convergence of 
product states. 

Lemma 11.15. If a CSCA t is fractal, so is t" for n e N. 

Proof. We prove that i™ is fractal by showing that it can neither be periodic nor have gliders. Ob- 
viously, no power of t can be periodic if t is not periodic, so only the glider case remains. If t n has 
minimal gliders £, Z then by 

t n (tz) = tt n z - X(tZ) 



12 



Johannes Giitschow, Sonja Uphoff, Reinhard F. Werner, and Zoltan Zimboras 



t£ is also a glider for t n with the same eigenvalue A and thus a multiple of £. Hence t£ = A£ holds 
for a monomial 3 A. If A = 1 the automaton would be periodic by Proposition III. 14l which is already 
ruled out. For A = w ±n , which is the only other possibility, t has gliders. Thus for t n to have gliders 
t has to have gliders. This is a contradiction to the assumption that t is fractal. So any power of a 
fractal CSCA is always fractal. □ 

III. INVARIANT STATES AND CONVERGENCE 

In this section we will consider different types of states on the spin chain and their evolution 
under CQCA action. Our focus lies on the search for invariant states and the convergence of other 
states towards invariant states. We consider special sets of states: product states, stabilizer states 
and quasifree states. Because CQCAs act in a translation-invariant manner, it is natural to look at 
translation-invariant states. We will therefore only consider those. 

The only state that we know to be invariant for all CQCAs (and all QCAs) is the tracial state, 
which vanishes on all finite Pauli products except the identity and is defined as the limit of states that 
have the density operator |l for each tensor factor. We strongly suspect this state to be the only 
invariant state for fractal CQCAs, but have no complete proof yet. 

A. Invariant states for periodic CQCAs 

It is obvious, that all states are periodic under the action of periodic CQCAs. Therefore a state 
is cither invariant or does not converge at all. Finding invariant states for periodic automata is in 
general very easy. For example, if a CQCA satisfies T p = 1 for some finite p € IN + , then all states 
of the type ^ Xm=i w ° are invariant. Moreover, any invariant state is of this form, because for uj 
invariant w = i Y^n=i w ° '• Finding pure invariant states is more complicated. In Section IlII C II 
we show, that for some period-two CQCAs pure invariant stabilizer states exist. For CQCAs without 
propagation there also exist invariant product states. Namely, if one Pauli matrix is left invariant by 
such a CQCA, then the state that gives expectation value one on this matrix and vanishes on the 
others is left invariant by this CQCA. A state with the same expectation value for all Pauli matrices 
is always left invariant by a non-propagating CQCA. 

B. Invariance and convergence of product states for non-periodic automata 

In this section we will only consider states which are translation-invariant product states with 
respect to single cell systems. They are of the form w(w(£)) = Yii 3K w (£i))> where $ is a state on a 
single cell. 

Proposition III.l. For a glider CQCA T there exist no translation-invariant product states that are 
T-invariant except the tracial state. All other translation-invariant product states weakly converge to 
T -invariant states 

Woo = limwo T". (18) 

n— too 

with the following properties: 

. w oo (w(0) = l ) ffw(0 = l, 

• |o; 00 (w(^))| < 1, &/w(£) is a product of gliders, 

• Woo(w(£)) = otherwise. 

If = 1 for some j, then uioo is the tracial state. 



3 By Remark III. 9 1 the only possible eigenvalues are monomials. 
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Proof. First we prove the non-existence of invariant product states. A state is invariant if u> o T = u>, 
i.e., u>(TA) — uj(A), VA g 21 has to be fulfilled. We require the automaton to be non-periodic which 
implies a finite (non-zero) propagation, so at least two 4 of the Pauli matrices are mapped to tensor 
products of at least three Pauli matrices 5 . Thus for these two Pauli matrices w(£) 

$(w(0) = w(w(0) 

= w(rw(0) 

holds. For all Pauli matrices <7j, |$(<Tj)| < 1. If |$(<7i)| = 1, then $(<7j) = 0, Vj ^ i. Now lets assume 
that $(<7j) 7^ for some i. The image of a single cell observable has to include at least two different 
types of Pauli matrices (different from the identity, e.g. o\ and 172 )• Else £+ = or £_ = or £ + = 
each case implying common divisors. 

Let us consider the case when there exists a Pauli matrix which is not mapped to a tensor product. 
It can not be mapped to itself, because then by Proposition lIL14l the automaton would be periodic. So 
it has to be mapped to another Pauli matrix which has to expand in the next step. Thus we only need 
to consider the case of expanding Pauli matrices. The image has to consist of more than one kind of 
Pauli matrices, so |$(erj)| = 1 is ruled out as an invariant state. Moreover, the image can not contain 
the original Pauli matrix even once, because 3>(<7i fe ) = (IL J - e jV'\fc<I> {a ij ))$(crj s . ) implies |<J>(<t^)| = 1, Vij 
which is already ruled out. So no Pauli matrix may occur in its own image, particularly not in the 
central position. If we only consider those central positions, we get a local automaton. It is an easy 
calculation to show, that all of these automata, which map no Pauli matrix onto itself have trace 
one. So the trace of our CQCA contains a constant which is a contradiction to the condition that it 
has gliders. So $(fTi) = 0, i = 1,2,3 and the only invariant state is the one with the density matrix 
p$ = il, i.e., the tracial state. 

In the following we consider the convergence properties of product states. Clearly, it suffices to 
establish the convergence of Pauli products. First, suppose that some |<E>(erj)| = 1. Then, according 
to Lemma [A. 11 for all Pauli products different from 1, and all times except at most one, the evolved 
product will contain a Pauli matrix different from o~j, and hence will have zero expectation in uj o T n . 
Hence uj^ is the tracial state. 

To treat the remaining cases, we hence assume from now on that |$(c7j)| < A < 1 for all three j. 
Hence if some Pauli product has k factors different from 1 , its expectation in oj is at most A fc . So let 
k(t) be the number of non-identity factors in the t th iterate of some Pauli product. If k(t) diverges, 
we have nothing to prove. So we may assume from now on that there is a constant k max < 00 such 
that k(t) < k max infinitely often. We focus on the subsequence for which this is the case. 

If the overall length of the Pauli product (largest-smallest degree) remains finite, and since we have 
assumed the absence of periodic finite configurations, then we must have a glider for some power of 
T. As argued in the proof of Lemma Til. 151 this is also a glider for T. In fact, for any product of glider 
elements, the left going and the right going gliders will eventually be separated, and from that point 
onwards the w-expectation will not change anymore. Hence the limit does exist and will be some 
number of modulus < 1. 

Hence we need only consider the case that k max is finite, but the positions of non-identity Pauli 
factors get more and more spread out. This is only possible, if some Pauli products near the edges 
of the given product have a similar property: no sub-product which gets widely separated from the 
rest is allowed to develop configurations with unbounded k max , because then the overall bound could 
not hold. Thus we again find bounded configurations with smaller k max , which may again be smaller 
gliders, or split up even further. By downwards induction we just get a complete decomposition into 
gliders for any configuration with non-divergent k rnax . This completes our proof. □ 

Proposition III. 2. Under the action of fractal CQCAs all product states converge to the tracial state. 

Proof. To show convergence for the fractal case, we need the results from Section III B 31 and Ap- 
pendix IA 11 They state, that a tensor product with only one kind of Pauli matrices can occur only 



4 The image of the third Pauli matrix is always determined by the product of the other two. 

5 The image has to be a tensor product of at least tree Pauli matrices, because an identity in the middle is not allowed. 
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once in the history of a fractal CQCA and that the number of non-identity tensor factors grows un- 
bounded. With the arguments used in the proof of Theorem III. 1 this means that for fractal CQCAs 
any given product state converges to the state which gives zero expectation value for any non-trivial 
tensor products of Pauli operators, i.e., to the tracial state. □ 



C. Invariance and convergence of stabilizer states 

In this chapter we consider pure translation-invariant stabilizer states. A stabilizer state is the 
common eigenstate of an abelian group of operators (usually tensor products of Pauli matrices) called 
the stabilizer group S. For those not familiar with the stabilizer formalism we give a short example. 

Example III. 3. The stabilizer group S — {a\®a\, 03&X73) stabilizes the Bell state ip = l/\/2(|l, 1) + 
|0, 0)). We check this by applying the stabilizer generators to the state (normalization is omitted): 

(Ol ® Ol)V> = (<7i®(Ti)(|l)®|l)) + (CTl®(Tl)(|0)(g)|0)) 

= |o,o) + |i,i) = v 

(<73®0- 3 )V = (cr 3 CT 3 )(|1) (g) |1)) + (ff 3 ® <7 3 )(|0) ® |0)) 

= (_|i))®(-|i)) + |o)®|0> 

= (-i) 2 |i,i) + |o,o) = v 



Now we move on to translation-invariant stabilizer states. Consider the stabilizer group S = 
(w(-P£), x € Z) generated by the translates of a Pauli tensor product w(£), where t x denotes the x-th. 
power of the phase space translation. The unique state u satisfying uj(A) — 1 for all A € S is called 
the pure translation-invariant stabilizer state corresponding to the stabilizer S Q. In the following 
we will refer to the stabilizer group simply as "stabilizer" . 



1. Invariance 

Let us first search which CQCAs leave pure translation-invariant stabilizer states invariant. 

Proposition III. 4. Only periodic CQCAs can leave pure translation-invariant stabilizer states in- 
variant. For each such stabilizer state with stabilizer S = (w(f x £), x £ Z) the CQCAs that leave the 
state invariant are periodic and form a group. The corresponding CSCAs are 



l + < 2 
a£ 1 + 



for a centered palindrome a. 

Proof. The invariance condition for a stabilizer state with stabilizer S = {w(f x ^), x G Z) is 

i£ = a£,a€ V. 

In i it was proved that for every translation-invariant stabilizer state u with stabilizer S = 
(w(f x ^), x <G Z) there exists a CQCA B, which maps w(0, 1) to w(£). Thus the condition becomes 

tb©=«b©. 

If we know automata that leave the "all spins up" state invariant, we can construct automata for 
arbitrary translation-invariant pure stabilizer states via 

*© = <) 
& btb-^Q = ab(°) 
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The only type of CQCAs that leave the "all spins up" state invariant are the shear transformations, 
which are represented by matrices 




with some palindrome a. A direct computation shows that the factor a has to be invertible. As 
we are in characteristic two and deal with centered automata, only a = 1 is possible. For a general 
translation-invariant stabilizer state the CQCAs that leave it invariant are represented by the CSCAs 
with matrices 



t c (o) 



i + 




The product of two such CQCAs also leaves the state invariant with t^(a)t^(b) — if (a + b). The 
inverse of each periodic CQCA is one of its powers, so it is also included in this set. We therefore 
have a group of period two CQCAs for each translation-invariant pure stabilizer state that leave this 
state invariant. □ 



2. Convergence 

We now consider an arbitrary translation-invariant stabilizer state w w (£) and the respective stabilizer 
S = (w(f a: £), x e Z). There exists a CQCA B satisfying B(w(£)) = 03. We rewrite our state in 
terms of the (73-state 

^w(c) = ^t 3 B^ 1 (2°) 

and use 

= w n (B- 1 T n ($$<r i )) 

= u n {{B- 1 TB) n B- i (®a i )). 

B~ 1 TB is also a CQCA of the same type (same trace) as T and B~ 1 (($Z) tr») is a tensor product of 
Pauli matrices. uj a3 is a product state and thus converges according to Propositions IIII. 1 1 and IIII.21 for 
glider and fractal CQCAs. 



D. Stationary quasifree states and convergence of quasifree states for one-step glider CQCAs 

In the previous sections we discussed stationary states and convergence of states under general 
CQCA actions. In this section we consider the particular glider CQCA G which is represented by the 
CSCA 

\1 It + 1 + u I 

because for this CQCA we can obtain new types of stationary states and new convergence results by 
employing the Araki- Jordan- Wigner transformation. Using the Theorem for glider equivalence III. 121 
we can construct invariant states for all automata with gliders that move one step in space every 
timestep. According to this theorem for any speed-one glider CQCA B, there is a CQCA A such that 
B = AGA^ 1 , and if is a G-invariant state, then u> o A^ 1 will be £?-invariant. 
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1. Araki-Jordan-Wigner transformation 

The Jordan- Wigner transformation is a way to map a finite spin-chain algebra A4f N to the algebra 
of a finite fermion chain. This method has been extensively used in solid state physics [liM2ll |. 
However, the method cannot be carried over directly to two-sided infinite chains. One has to introduce 
an additional infinite "tail-element" for the transformation to work. This extended transformation 
was introduced by Araki in his study of the two-sided infinite XY-chain [2^ |. and it is sometimes 
referred to as the Araki-Jordan-Wigner construction. 

The C*-algebra describing a two-sided infinite fermion chain is the algebra J-"=CAR(^ 2 (Z)), i.e., it is 
the C*-algebra generated by 1 and and the annihilation and creation operators {c x } xe z and {c*} xe z, 
satisfying the canonical anticommutation relations: 

C x C y + C y C x ~ ^x,s/l) C x C y + C y C x = ®- 

The translation automorphism rjr on this algebra is defined by Tjr(c x ) — c x+ \. J- is isomorphic to 
the observable algebra 21 of the spin chain, but there exist no isomorphism i: 21 — > T that satisfies 
the property i o r = rjr o i. This intertwining property would be needed to derive the translation 
invariance of a state u o i on 21 from that of 10 on T . This problem can be circumvented by the 
Araki-Jordan-Wigner construction. 

The ordinary Jordan- Wigner isomorphism between the iV-site spin-chain algebra M® N (generated 
by the finite number of Pauli matrices {af , erf , erf }xe{o,i,...,N-i} ) 5 and the iV-site fermion-chain 
algebra (generated by {c x , c*} xe { ,i,...,N-i}), is given by 

x-l 

LjwK) ■= \{{2c y c* y -l){c* x + c x ), 

y=0 



x-1 



,N I „x 



(trf) := U(2c yC ; - 



L jw( a 3) '■— 2c x c x 



However, as we have mentioned, the Jordan- Wigner transformation cannot be generalized to be a 
translation-intertwining isomorphism between the two-sided infinite spin and fermion chains. In an 
informal way, one could say that an element of the form "IIy=-oo(^ c y c y — wom d be needed in the 
definition of a "two-sided infinite chain Jordan- Wigner transformation" . However, T doesn't contain 
such an element. The basic idea of the Araki-Jordan-Wigner construction [22] is to extend the algebra 
T with such an infinite tail-element. More concretely, one defines the C*-algebra F to be the extension 
of J 7 by a self-adjoint unitary element U satisfying: 6 



Uc x U 



c x if x > 
—c x if x < 



Clearly, every element of J- can be uniquely written in the form a + Ub with a, b E J 7 , i.e. T = F+UJ-. 
One can extend the translation automorphism tjt to J- through the formula fjr(a + Ub) :— Tjr(a) + 
C/(2coCq — V)Tjr{b). Let us define the elements 

(UyZl^Cyc; - 1) if ^ > o 

A x :=ll if x = , 

Ini^vs-i) if ^ < o 

and let us introduce the elements 

c?T ~UA x (c*+c x ), o? :=UA x i(c* x -c x ), ^ := (2c x c* - 1). (22) 



In the literature the symbol T is used almost exclusively for denoting this unitary element. However, we chose to 
denote it by U to avoid confusion with the time-evolution automorphism, which is denoted by T in this paper. 
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Denote by T the C*-subalgebra of T which is generated by the elements of {a\ x , , a^ x \ x £ Z}. A 
direct computation shows that the elements defined in (|2"2"j) having the same "spatial index" x satisfy 
the Pauli-relations, while any two of these elements having two different spatial indices commute. 
Hence J- and 21 are isomorphic, and an isomorphism is given by the map II : J 7 — > 21, defined as 

n(^T) = a x , n(SJ x ) = cxf , n(^ x ) = <rf , Vx e z. 
Moreover, if we denote by fjr the restriction of fjr to T then 

T O II = II O Tjr, 

i.e., II intertwines the translations of the two algebras. 

Let w be a translation-invariant state on the fermion-chain algebra J- ' . By defining uj(a+Ub) := 10(a) 
we get a translation-invariant state on J-. Restricting this state to J- we get a Tjr-invariant state cD, 
and uj jw =Q oil 1 will be a translation-invariant state on the quantum spin-chain. In this way we 
can transfer translation-invariant states from the fermion-chain to the spin-chain. 

Any CQCA automorphism T can naturally be transferred to an automorphism on T commuting 
with Tjr by the definition T := no To II -1 . In the case of the glider CQCA G (|2"Tj) we can do even 
more. The transferred automorphism G (characterized by G{a\ x ) = o^ x , G(a 2 x ) = erj 1 1 ~0\ X 02 X+1 ) 
can also be extended to an automorphism G : F — > F in a translation- invariant way (Gof F = f F oG) 
with the following definition: 8 

G(U) = iU(c_ 1 + c*_ 1 )(c -4). 

It is exactly this type of "translation-invariant extension property" that allows us to find stationary 
translation-invariant states of the glider CQCA by the Araki- Jordan- Wigner method. 

Restricting the G automorphism to the fermion-chain subalgebra J- C J- , we obtain the automor- 
phism Gjr, which acts in the following way: 

Gjr(c x ) = -(c*_ x + c x _ x + c* x+1 - c x+1 ). 

The automorphism Gjr : J- — > J- takes an especially simple form in terms of majorana operators. 
These operators are defined as 

m 2x := i(c x - c*), m 2x+1 := c x + c*, (23) 

for any and they generate F. The action of Gjr on these operators is 

Gjr(m 2x ) = m 2x -2, Gjr(m 2x+ i) = m 2x+3 . 

Clearly, if we find a state w on J that is both Tjr - and G^-invariant, then the Araki- Jordan- 
Wigner transformed state uj jw will be a r- and G-invariant state on the quantum spin-chain. In 
the next section, we will recall the definition of quasifree states on fermion-chains and then deter- 
mine the translation- and Gjr-invariant quasifree states. In this way, using the Araki- Jordan- Wigner 
construction, we can obtain a whole class of translation- and G-invariant states on the spin-chain. 

2. Stationary quasifree states 

A state u> : J- —¥ C is called quasifree if it vanishes on odd monomials of majorana operators 

ui(m xi . ■■m X2n _ 1 ) = 0, 



7 It is clear that by this definition u will be a normalized functional on F. The positivity of ui follows from ui({a + 
Ub)*{a + Ub)) =aj(a*a + 6*&) > 0. 

8 We only have to define the image of U under G, since any element / £ T can be uniquely written as a linear 
combination / = f% + Uf2 , where /i , G T . 
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while on even monomials of majorana operators it factorizes in the following form: 

n 

uj(m Xl ...m X2n ) =^sgn(Tr)Y[u(rn Xni2l _ 1) m Xnm ), 

7T 1 = 1 

where the sum runs over all pairings of the set {1,2,..., 2n}, i.e., over all the it permutations of the 
2n elements which satisfy 7r(2Z — 1) < ir(2l) and tt(21 — 1) < n(2l + 1). Hence, if we assume that 
Xi 7^ Xj when i ^ j, then ui{m Xl . . . ni X27l ) is simply the Pfaffian of the 2n x 2n antisymmetric matrix 
Aij := Uj(m Xi m Xj ). 

An automorphism a : J- — > J- that maps any majorana operator onto a linear combination of ma- 
jorana operators is called a Bogoliubov automorphism [23| . For any quasifree state uj, the Bogoliubov 
transformed state uj' :— uj o a will again be quasifree. 

A quasifree state uj is translation-invariant, i.e. uj o tj- = uj, iff uj(m x m y ) — uj(m x+2 m y+2 ) for all 
x, y G Z. Translation-invariant quasifree states are characterized by a majorana two-point matrix 
which is a 2 x 2-block Toeplitz matrix of the form (for a proof, see e.g. [24| ) 



uj(m 2x m 2y ) U}(m 2x m2y + i) 
u{m 2x+1 m, 2y ) u)(m 2x+1 m 2y+1 ) 



1 

2^ 



-*p(*-i/)dp, 



(24) 



where are L°°([— 7r, it]) functions, and the matrix function 



eL|= x2 ([-7r,7r]) 



satisfies 



(Q^&O^Q^Cp). g (w) (-p) = 2-l-(Q^(p)) T , 0<Q^(p)<2.1 



(25) 



almost everywhere (here (Q^(p)) T denotes the transpose of Q^(p)). A translation-invariant quasi- 
free state w is pure iff for almost every p the eigenvalues of Q^\p) are either or 2. Q^(p) is called 
the symbol of the majorana two-point matrix = uj(m x m v ). 

Now we are ready to characterize the translation-invariant quasifree states that are stationary with 
respect to the time-evolution Gjr (defined in the previous subsection). 

Proposition III. 5. A translation-invariant quasifree state uj is invariant under the Gjr automorphism 
iff the symbol of its majorana two-point matrix has the following form: 



q$(p) 



where and q^ 2 are real L°°{[— n, tt]) functions that take values between and 2 (almost every- 
where). 

Proof. The Gjr automorphism acts on the majorana fermions as 

Gjr(m 2x ) = m 2x -2 , Gjr(m 2x+1 ) = m 2x+3 , 

hence it is a Bogoliubov automorphism. Moreover, Gjr commutes with the translations. Thus uj' — 
uj o Gjr will again be a translation-invariant quasifree state, and uj' is equal to uj iff its majorana 
two-point matrix is the same as that of uj. The majorana two-point function of uj' is 



lo' (m 2x m 2y ) uj' (m 2x m 2y+1 ) 
lo' (m 2x+1 m 2y ) lo' (m 2x+1 m 2y+1 ) 



uj(Gj-(m 2x m 2y )) uj(Gj-(m 2x m 2y+ i)) 
u}{Gj-{m 2x+1 m 2y )) u}(Gj-{m 2x+1 m 2y+1 )) 

u(m 2x ^ 2 ni 2 y- 2 ) uj(m 2x - 2 m 2y+3 ) 
u{m 2x + 3 m 2y - 2 ) u(m 2x+3 m 2y+3 ,) 

uj(m 2x m 2y ) uj(m 2x - 2 m 2y+3 ) 
u(m 2x+3 m 2y ^ 2 ) uj(m 2x+ im 2y+1 ) 
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where we have used that uj is Tjr-invariant, and that Tjr(m x ) = m x +2, which follows from the definition 
of the majorana operators ()23[) . Comparing the majorana two-point functions, one can see that 
uj' = w iff uj(jn2xm2y+i) = u{m2 X -2'm2y+3) for any x, y € Z. From form (|24j) of the majorana two- 



point functions of translation-invariant quasifree states we know that there exists a q^\ 
function such that 



-7T, 7T\ 



IT 

u(m 2x m 2y+1 ) = i- J q { $(p)e-^-vUp. 

— 7T 

The condition oj(m2x'^2y+i) = ^(^2^-2^2^+3), which should hold for all x,y £ Z, means that 

7T 7T 



2tt 



9i2 (p)(1 - e^e-^-^dp = 



must be satisfied. It is a well known theorem in functional analysis, that an L°°([— it, it]) function for 
which the Fourier transformation vanishes must be almost everywhere zero, hence q{^(p)-(l-e 2ip ) = 0, 



from which one concludes that q^ip) = 0- This means that for a tjt- and Gjr-invariant state uj the 
symbol of the two-point majorana matrix has to be of the form 



«i2(p) 



where q^ and q^ij are L°°([— 7T, 7r]) functions, and according to ([23)1 they also have to satisfy the 



» 



inequalities < gjy (p) < 2 and < q^ (p) < 2 almost everywhere. Thus we have arrived at our 
proposition. □ 

Using the Araki- Jordan- Wigner transformation we can transfer such a tjt- and Gjr-invariant qua- 
sifree state u> on T to a state uj jw on the spin chain which is r-invariant and stationary with respect 
to the glider CQCA G. 



3. Convergence of quasifree states 

In this section we will show that under the repeated action of the Gjr automorphism, any translation- 
invariant quasifree state will converge to one of the Gp- invariant states specified in Proposition MI. 5l 

Proposition III. 6. Let uj : J- — > C be a translation-invariant quasifree state with a majorana two- 
point matrix that belongs to the symbol 



? M (p) = 



&p) &P) 



where Q (w) eL^ 2 ([ — 7T,7r]) satisfies the relations \25\) . The n time-step evolved state is denoted by 
uj n := uj o G 7 ^. The pointwise limit uj^A) := lhxin^oo uj n (A) exists for all A G T, and the function 
Woo : J~ —t C defined in this way will be a translation-invariant quasifree state with a majorana 
two-point matrix that belongs to the following symbol 



Q (w ~>(p) = 



g£J(p) o 

<$/(?) 
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Proof. In the proof of Proposition IIII.51 we showed that if is a translation- invariant quasifree state, 
then W2 = wo Gjr will also be such a state. By induction it follows that w„ is also a translation- 
invariant quasifree state for any n E N + . Hence for an arbitrary odd monomial of majorana operators 
u> n (m Xl m X2 ■ ■ -m X2 

jv+i) — Oj and the limit ui(m Xl m X2 ■ ■ ■ ni X2N+1 ) :— lim„_ ! . 00 co n (m Xl m X2 ■ ■ ■ itl X2N+1 ) 

exists and is zero. 

Next, we prove the pointwise convergence of the majorana two-point matrix. The majorana two- 
point matrix of u n is 



Lo n (m 2x m 2 y) u n {m 2x m 2y+ i) 
u n (m 2x+ im 2y ) Lo n {m 2x+ im 2y+1 ) 



u(G n T {m 2x m 2y )) uj(G^(m 2x m 2y+ i)) 
u(G n jr{m 2x+1 m 2y )) u{G^{m 2x+1 m 2y+1 )) 

u(yn 2x - 2n rn 2 y- 2n ) (jj(m 2x - 2n m 2y+ i +2n ) 
w(m 2x +i+ 2n m 2y - 2n ) Lo{m 2x+ i +2n m 2y+ i +2n ) 

Lo{m 2x m 2y ) u(m 2x - 2n m 2y+ i +2n ) 
u{m 2x+ i +2n m 2y ^ 2n ) uj(m 2x+ im 2y+ i) 



this means that the symbol Q( Wn > is given by 



(?) 



The limits of elements of the two-point majorana matrix are: 



^oo{jn 2x m 2y ) Uoo{m 2x m 2y+ i) 
Uoo{m 2x+ im 2y ) w 00 (m2s+im2y+i) 



1 

27 

1 

2^ 



lim uj n {m 2x m 2y ) lim Lo n (m 2x m 2y+ i) 

n— >oo n— ¥oo 

lim io n (m 2x+ im 2y ) lim Lo n (m 2x+ im 2y+ i) 

lim r q[^{p)e- l <- x -y^dp lim T q^ (p)e~ i ^y~ 2n ^dp 



lim C q 2 L °J(p)e- i( - x -y+ 2n ')Pdp lim T q\ U J Me-^^dp 



flj?(p) 
o 



«& } (p) 



-ip( x -y) 



dp, 



(26) 



where we have used the Riemann-Lebesgue lemma, which states that for any integrable function / 
defined on the interval [a, b]: 



lim / f(x)e lzx dx = 0. 

a 

For an arbitrary even monomial of majorana operators the convergence can be proved by: 
Uoo(m Xl m X2 ■ ■ ■ m X2 

at+i) — li m w n ijn Xl m X2 ■■■ m x 

2N + l) 

n— >oo 

JV 



= lim Vsgn(7r)TTw„(m^ (!W) m t „ (3i) ) = 

iV 

= sgn ( 7r ) n I ^ w ™( m »-(«-i) m ^(2!)) 

* — » A A Lre— >oo ^ ' * ' 

JV 



(27) 



Since the finite linear combinations of majorana operators form a norm-dense subset in T the limit 
lirn„_ i . 00 w„(/) must exist for any f € T due to the uniform boundedness of the states in the sequence, 
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and hence we have obtained with this pointwise limit a linear functional Woo : J- — > K. Moreover, 
ujoo is uniquely determined by the values it takes on monomials of majorana operators, and since 
the equations (|26j). ([27]) are satisfied can only be the quasifree state for which the symbol of the 
majorana two-point function is given by (|26[) . □ 

It is worth mentioning that although a pure quasifree state will, of course, stay pure for any time 
t, taking the weak, i.e. pointwise, limit of the states when t —¥ oo one obtains a mixed state (if the 
original state was not G- invariant). In quantum many body physics such relaxation from a pure to 
mixed states has been studied in a quite different setting, namely how certain time averages of pure 
states that evolve under a Hamiltonian dynamics can be described by mixed states (25|. It would be 
interesting to look at such weak limits when t — > oo also for Hamiltonian evolution (for systems with 
infinite degrees of freedom), and to look at the connection between these two approaches. 

Finally we have to point out, that the results of this section using the AJW transformation can 
generally not be transferred to other CQCAs than the glider CQCA G, because most CQCAs don't 
have the form of Bogoliubov transformations on the CAR- Algebra. Automata with neighborhoods 
larger than nearest neighbors can map creation and annihilation operators to products of these on 
the CAR-algebra and are thus in general not Bogoliubov transformations. (An obvious exception are 
powers of G.) All automata, that don't leave at least one Pauli matrix locally invariant don't allow for 
a tail element that is left invariant under the transformed CQCA. These automata are characterized 
by a constant in their trace polynomial. In conclusion, only nearest neighbor automata without a 
constant on the trace, i.e. some glider and period two automata and their powers, can be transferred 
to Bogoliubov transformations on the CAR-algebra. For more details see [17j . 



IV. ENTANGLEMENT GENERATION 



In this section we investigate the entanglement generation properties of CQCAs. First we de- 
rive general bounds for the entanglement generation of arbitrary QCAs on the spin-chain both in a 
translation-invariant and a non-translation-invariant setting. Then we investigate CQCAs acting on 
stabilizer and quasifree states. We find that the entanglement generation is linear in time, similarly 
to the case of Hamiltonian time evolutions [29| . 



A. General bounds on the entanglement generation of QCA 

In this section we derive general bounds on the evolution of the entanglement of a finite number 
of consecutive spins with the rest of the chain under the action of a localized automorphism (e.g. a 
QCA). We will only consider the case when the whole chain is in a pure state. In this case the proper 
measure of entanglement is given by the von Neumann entropy 

S = -Trp s log 2 Ps, 
where ps is the reduced density matrix of the finite segment of spins. 9 

1. The non-translation-invariant case 
Theorem IV. 1. Consider the observable algebra of an infinite chain of d-level systems 

21^ := (g) 2lf, a^SM* 

i— — oo 



Our bounds for the entropy generation hold also for non-pure states, although in this case the von Neumann entropy 
is not directly related to entanglement. 
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and an automorphism T : 21^ — ^ 2t( d '. Let us introduce the notation: 21^ m2 \ '■= ® 7 k=m 1 ^fc^ ■ 
Suppose that T satisfies the locality condition 

W\l k2] ) c agy , 

for some fixed integers k\, fc 2 , h, fa, with l\ < k\ and k 2 < ^2- 

Lei uj be a state on the spin-chain % <yd \ and let us define the T -evolved state as oj' := w o T. 

TTie restrictions of these states to a subalgebra 21^ m ^ will be denoted by W[ mi)m2 ] and wj mi m 1, 
respectively. Then the following bounds hold for the von Neumann entropies of the restricted states: 

S^[ kl ,k 2 ]) ~ 2n ^g 2 d < S(uj' [ki k2] ) <S(ui [klM] ) + 2n\og 2 d, (28) 

where n = fa — l\ — k 2 + k\ . Moreover, these bounds are sharp. 

Proof. Restricting the automorphism T to the subsubalgebra 21^ fc i, we get a monomorphism 10 
T[k 1: k 2 } '■ 21 [it j k 2 ] ~^ ^Ri I2V This monomorphism can be extended to an automorphism T : 2li^ m —> 
^Ih h]' 11 us Produce the following state on 2l[; lj ; 2 ]: 

U'([h,l 2 ]) := W^^] o T. 

From this construction it immediately follows that fa])[ki,ka] = fe 2 i- Moreover, since 

^'(['17^2])) and wjjj ; 2 ] are connected by an automorphism their von Neumann entropies are equal: 
S(w'([M 2 ])) = S( W[h , 2] ). 

We will prove the bounds (|28|) using the subadditivity of the von Neumann entropy. The subchain 



2l^ 2 j can be divided as 21^^ = 21^ k ^ ®2l^ t , where 2l^? s t is isomorphic to the algebra of d n x d n 

matrices, hence the maximal entropy of a state defined on 21^,* is nlog 2 d. The triangle inequality 
and the subadditivity theorem give the following inequalities: 

S ( u [ki,k3]) ~ n\og 2 d < S(u [luh] ) < S(u [klM] ) + n\og 2 d 
S HkiM])- nl °S2d< 5(w'([li,/ 2 ])) <5(w[ fclffca] )+nlog 2 d 



Now, using that S(uj[i lt i 2 ]) = S(uj'([h, fa])) we immediately obtain the bounds ([28]). 

The sharpness of the inequalities follows if we consider a state on the spin-chain where the sites at 
2i are maximally entangled with the sites at 2i + 1 and we consider the translation r which just shifts 
all one-cell algebras by one cell to the right as our time-evolution. Then k 2 — ki = 2 and fa — h = 3, 
and we get n = fa — fa — k 2 + fci = 1. Now, restricting this state to the subalgebra 2l 2 ^ 2i , 3 the 
entropy of the restriction is zero. However, the entropy of this restriction after the time evolution 
will be 2 log 2 d, since the two sites at the border will be maximally entangled with sites outside the 
considered region. Let us note, that in some sense we broke the translation-invariance only minimally, 
since the considered state is invariant under the square of the translations r 2 . 

In the above example the generated entanglement is destroyed in the next step, so the bound is 
only saturated for one time step. But a slightly more involved example shows, that the bound can be 
saturated for arbitrarily many timesteps: 

Let again T be the translation automorphism on 2l( d * 1 , and let us consider the state on a spin-chain 
which is defined as the direct product of totally mixed states between the lattice site at i and the 
lattice site at — i + 1 for all i. So in this state the lattice site at 1 is fully entangled with the lattice site 
at 0, the lattice site at 2 is fully entangled with the lattice site at —1, and so on. Now, we will consider 
subsystems of arbitrary length k = k 2 — k\. If k is even, k — 2j, then the subsystem we consider is the 
interval [—j + Its original entropy is 0, and the entropy grows linearly during the time-evolution 
saturating our linear bound until it reaches the maximal entropy it can obtain, namely k log 2 d. After 
this it stays constant. If k is odd, k = 2j + 1, we consider the interval as our subsystem. The 

original entropy of the subsystem is log 2 d, and the entropy grows linearly and saturates our bound 
until it reaches (2j + 1) log 2 d, then it stays constant. □ 



An injective but not (necessarily) surjective homomorphism. 

This can be simply seen by noting that h] = T(9t[^ ® OS, where <B ((g^LT, 1 ^) ® «8>j 2 =fc 2 + i a j d> )> if Q 
is an isomorphism between the latter two algebras, then one can define Tn lt i 2 ] as 7] fcl fe2 ] <g) Q. 
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2. The translation-invariant case 



In the previous subsection it was shown that the bounds (|28[) on entanglement generation are sharp 
in the general case. However, considering translation- invariant states and QCA automorphism, i.e., 
automorphisms that commute with the translations, we can sharpen these bounds further. 

Theorem IV. 2. Consider the observable algebra of an infinite chain of d-level systems 

+00 



2t (d) := (g) 2lf, 2lf-A/ d 



and a QCA automorphism T acting on having a neighborhood of n "extra cells", i.e., T is an 
automorphism satisfying 

T(2l,) (d) C 2l^ ni ® 2t^ ni+1 ® • • • ® Sag? na _ 1 ® ag^ , T o t = t o T , (29) 

where i eZ, t is i/ie translation automorphism ontyS d > , andn\ andn-i are integers satisfying ni+n-2 = 
n>0. 

Let oj be a translation- invariant state on the spin-chain, and let us define the time-evolved state 
(at time t € N) as oj(t) := uj o T . TTie won Neumann entropy Sx,(t) 0/ i/ie restriction of tu(t) to L 
consecutive qubits can be bounded in the following way: 

S L (0) - nt log 2 d < 5 L (i) < (0) + rrf log 2 d. (30) 

Moreover, these bounds are sharp for d = 2. 

Proof. Since the state u is translation invariant and T commutes with the translations, the "entropy 
production" is the same for the automorphisms T and Tot' {k 6 Z), hence we can assume without 
loss of generality that in Eq. (|2"9"|) m, 712 > 0. 

Let us denote the restriction of a state (p : 21^ — > C to 21^ , = ®£!? mi 21^ by ^j^j, 
Consider the subalgebra 2l[oL_i] °f 2l^ d \ which corresponds to L qubits. Restricting the auto- 
morphism T to this subalgebra, we get a monomorphism : Stf^i— 1] ~ ^ ^[-ni i-i+n 3 ]' ^ e 
will also consider the inverse automorphism T _1 , and restrict to a monomorphism (T _1 )l : 
21[ol-i] — ^[-n 2 i-i+ni]' 12 ' ^^ ie monomOT pld sm can be extended to an automorphism 

f-l . «(<*) 9 ,(d) 

■ <*-[_ na ,X,_l+ ni ] ^ ^[-na.L-l+m]- 

Let w be a translation-invariant state on 21^, then cj(1) :=woT will be translation-invariant, too. 
Let us also define the following state on 21^ r , . , 



w([-n 2 , L-1 + rai]) := w(l)[_ tlaii _ 1+ni ] o T 



-1 

L ' 



The von Neumann entropy of uj([—n2,L — 1 + m]) and u>(l)i n2 x-i+nx] are t ne same (since they 
are connected by an automorphism), and it follows from the definition of T^ 1 that ri2,L — 1 + 
"•i])[o,i-i] = W [0,L-1]- 

Now, from the strong subadditivity of the von Neumann entropy it follows that for a translation- 
invariant state oj S(w[ mi . m2 j) > S(ui[k 1 ,k 2 ]) if 771 2 — m i > &2 — &i |28|, hence 

S L (1) = S(w(1) [0 ,l-i]) < S(w(l) [ _ nai£ _ 1+ni] ). (31) 

On the other hand, using the subadditivity of the entropy for the state u3([— 112, L— 1+ni]) (dividing the 

observable algebra of the subchain [— n 2 , ni] as: 2tf_ na L-i+m] = ^[0 l-i] ® (^[-n 2 -i]®^[i i-i+ni]))' 
we get: 

S L (u([-n 2 ,L -1 + m})) < SrL(w([-n2,L-l + n 1 ])[o,i-i]) + (ni + n 2 )Iog 2 d = 5i(0) + nIog 2 d (32) 



The fact that the range of 21[o,l-1] under the the action of T 1 is in 2l[_ n2 ,L—i+ni] was shown in [27jl . 
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Combing the fact that S , (w(l)r_ n2i £_ 1+ni ]) = S(cj([— ri2,L — 1 +ni])) with the inequalities (j3Tj) and 
(|32[) we arrive at the S(l) < S(0) + n log 2 d. By simple induction we obtain the desired upper bound: 

S L {t) < S L (0)+nt\og 2 d. 

The lower bound in (|30[) can simply be obtained by "reversing the time arrow": suppose that for a 
QCA automorphism T this lower bound does not hold, this would mean that for the QCA T _1 the 
upper bound wouldn't hold, which is a contradiction as we proved the upper bound just now. 

The sharpness of the inequalities for d = 2 follows from the study of of Clifford QCAs acting on the 
"all spin up state" in Section HVBI □ 



B. Entanglement generation starting from translation-invariant stabilizer states 

In this section we will consider the entanglement generation of CQCA acting on translation-invariant 
pure stabilizer states. We will first calculate the bipartite entanglement in a general translation- 
invariant pure stabilizer state. Using this result we will present a proof of asymptotically linear 
growth of entanglement for non-periodic CQCAs. 

For every translation-invariant stabilizer state uj with stabilizer S = (w(t x £), i£Z) there exists 
a CQCA T which maps w(0, 1) to w(£). Each S = (w(t x £), x E Z) defines a unique translation- 
invariant stabilizer state if and only if £ is reflection invariant and gcd(£+,£_) = 1 9|. The image 
of a one-site Pauli matrix under the action of a CQCA B is always of this form. The study of the 
entanglement generation of CQCAs acting on initially unentangled stabilizer product states is thus 
equivalent to the study of the entanglement properties of translation-invariant stabilizer states. 

There are several results on the entanglement entropy for stabilizer states in the literature, the most 
general example would be the formalism developed in [29J . One case considered is a bipartite split of 
the state loab with respect to the subsystems A and B. The set of stabilizers S then splits up into 
three sets. Sa and Sb are the local stabilizers, which act non-trivially only on part A resp. B. The 
third set Sab accounts for correlations between the subsystems. It is defined as follows: 

Definition IV. 3. The correlation subgroup Sab for a bipartite stabilizer state is generated by all 
stabilizer generators that have support on both parts of the system. 

Sa and Sb together form the so called local subgroup. The correlation subgroup Sab can be brought 
into a form where it consists of pairs of stabilizers whose projections on A (and B) anticommute, but 
commute with all elements of other pairs and the local subgroup. The entanglement or von Neumann 
entropy of such a stabilizer state is 

E{oj A b) = l\S AB \, (33) 

if loab is a pure state, where |<Sab| is the size of the minimal generating set of Sab ■ 

Unfortunately, in [29] only finitely many qubits are considered. The proof of (|33[) relies heavily 
on this fact. However, there is a different approach to the bipartite entanglement in stabilizer states 
which we use to extend this result to infinitely many qubits. 

In our approach we make use of the phase space description of stabilizer states introduced in Q . A 
stabilizer state is fixed by a set of defining commuting Pauli products. In the phase space description 
commutation relations are encoded in the symplectic form cr(£, if). If [w(£), w(ry)] = we have 
cr(£, 77) = 0. Thus abelian algebras of Weyl operators (Pauli products) correspond to subspaces on 
which the symplectic form vanishes. Those subspaces are called isotropic subspaces. If for an isotropic 
subspace T o~(r),£) = 0, V( € I implies r\ G I, we call X maximally isotropic. Maximally isotropic 
subspaces correspond to maximally abelian algebras. In Q it was shown that the above condition on 
£ (reflection invariance and gcd(£_|_,£_) = 1) is equivalent to the condition that "P£ is a maximally 
isotropic subspace. V£, denotes the space generated by the products of £ and all elements of V . We 
have (w(r)),ri € V£) = (v?(t x £), x € Z) = S. VS, is the phase space of the stabilizer group S. By 
(~P£)a etc we denote the phase spaces of Sa etc. 

Theorem IV. 4. The number of maximally entangled qubit pairs in a pure translation invariant 
stabilizer state stabilized by S = (w(f :E £), x € Z) on a bipartite spin-chain is the number of pairs 

V 1 , C € (VOab untha A (v\C 3 ) = Sij, o-{rf A ,rf A ) = a(C,C j ) = 0, a{rf A ,n) = ^/Nm) = 0, V M e (VZ) A , 
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and a{C, A ,fjb) — cr(C l ,/i) = 0, V/i G (P£)a- where T)a, Qa denote the restriction of the phase space 
vectors to subsystem A completed with on B so we can use the the symplectic form a of the whole 
chain. 

Proof. If we restrict the stabilizer to subsystem A (or B) it is in general not translation invariant any 
more. Therefore the corresponding subspace is not maximally isotropic. The restricted state is not 
a pure translation invariant stabilizer state. However the uncut stabilizer operators in Sa stabilize a 
subspace of the statespace. The elements of the correlation subgroup Sab which is generated by the 
cut stabilizer generators map this subspace onto itself because they commute with the elements of Sa ■ 
From the theory of quantum error correction [30( we know that a pair of operators leaving a stabilized 
subspace invariant can be used to encode a logical qubit if the operators fulfill the same commutation 
relations as o~\ and 03. As the restrictions of the elements of Sab don't have to commute such pairs 
of operators can exist. In the phase space description this means, that we have to find 77, Q G (V£;)ab 
with o-{t]aXa) = 1, a(r) A ,fXA) = o-(i], (j,) = 0, V/i G {VCiA, and a(( A , Ha) = c(O^) = 0, V> G (P£,)a- 
Several such pairs encode several qubits. Of course the operators from different pairs have to commute. 
Thus the qubits are encoded by pairs of operators (a\, a\) = (w(?7 l ), w(£ 4 )) whose phase space vectors 
fulfill a A (v\C J ) = Sa, a(rj A ,rf A ) = a(C,C j ) = 0, <r{rf A ,n) = <r{r}\n) = 0, V/i G (P^a, and a(C A ,fx) = 
o-(C>m) = 0, V/i G (V£)a- As we have cj^aXa) + o-(t) B ,(b) = <r(r)A + Vb,Ca + (b) = (r(r), C) = 
0, V?7, C S we know that o~{t\b, (b) = &(va,Ca), thus the restrictions of our operators to system 
A and B fulfill the same commutation relations. We therefore have pairs of operators of the form 
(af ® a 1 , ® & B )- Each such a pair encodes a Bell pair as seen in Example lIII.31 Thus the number 
of maximally entangled qubit pairs is the number of such pairs of operators. □ 

We now show that the number of qubit pairs is 4 |«Sab j - As mentioned above, only Weyl operators 
w(£) fulfilling certain conditions can span the stabilizer S = (w(f x £), x G Z) of a pure state. On the 
level of tensor products of Pauli matrices the above conditions have three important consequences that 
stem from the requirement for £ + and £_ to have no common divisors and to be reflection invariant: 

1. The length of the product has to be odd, because palindromes of even length are always divisible 
by (1 + u). We will write I = 2n + 1. 

2. The central element of the product can not be the identity. Else £ has the divisor + u). 

3. At least two different types of elements (both different from the identity) have to occur (e.g. o\ 
and 1T2 )■ Else £ + = or = or £ + = each case implying common divisors. 

If we make a bipartite cut 13 in our system, 2n stabilizers will be affected. All other operators are 
localized on one side of the cut, only those with localization on both sides are elements of Sab- If we 
could find k pairs of anticommuting operators in the projections of Sab on the right (or left) halfchain 
there would be k pairs of maximally entangled qubits. In fact we can always find k = n such pairs 
and thus k = ^|<5>ab|- 

Definition IV. 5. The bipartite entanglement E(lj^) of a translation-invariant stabilizer state is 
the number of maximally entangled qubit pairs with respect to any bipartite cut. 

Theorem IV. 6. A pure translation-invariant stabilizer state of stabilizer generator length 2n+ 1 
entangles n qubit pairs maximally with respect to any bipartite cut. 

length(Q - 1 

£ k) = 2 ■ ( ' 

For the proof of this theorem we refer to Appendix IA 21 

Now it remains to show how the stabilizer generator length evolves under the action of CQCAs, 
and to deduce the asymptotic entanglement generation rate. 



Due to the translation invariance all possible cuts are equivalent. 
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Definition IV. 7. The asymptotic entanglement generation rate from stabilizer states for CQCAs is 
defined as 

ATT 2. 

— = lim -Efa,t), (35) 

l\t t— >oo t 

where E(uj^,t) is the bipartite entanglement at time t. 

We will now prove the following theorem: 

Theorem IV. 8. The asymptotic entanglement generation per step of a general centered CQCA is the 
highest exponent in its trace polynomial, dg(trt). 

For the proof we need the following lemma: 

Lemma IV. 9. The length 2n + l of the stabilizer generator of a stabilizer state grows asymptotically 
with 

An 1 

— = lim -n(t,0 =dg(trt) (36) 

At t->oo t 

for any centered CQCA T and any translation invariant pure stabilizer state wj. 

Proof. We know that CQCAs map pure stabilizer states to pure stabilizer states and stabilizer gen- 
erators to stabilizer generators [9(. Thus the image of a pure stabilizer state with stabilizer 
S = (w(-P£), x € Z) under T* is again a pure stabilizer state. We can write £ = bu) and 
lj^ = oj(o.i) ° B. Therefore we can write the evolved state as uj t t^ = lo^ o T* = U(o,i) ° T l o B. 
The length of the stabilizer generator is determined by the highest order of the stabilizer generators 
polynomials, dg(£ ). Namely the stabilizer generator is of length 2 • dg(£) + 1. So we have to calculate 

dg(t t o = dg(* t bO). 

An arbitrary product of CSCAs can be written as JliLi The series (a;)i<f<fc = dg(rj' =1 ti) is 
subadditive, i.e. a„+ m < a n +a m , because concatenation of CSCAs is essentially the multiplication 
and addition of polynomials which is subadditive in the exponents. For subadditive series a n Fekete's 
Lemma [3ljj states that the limit linin^oo ^ exists. In our case the series is always positive, so the 
limit is positive and finite. To determine the limit, we use the subsequence of the t = 2 k th steps. 
Using the Cayley-Hamilton theorem we get 

t 2 " = t(trtf- 1 +tY(trtf- 21 



i=l 



and 

:= dg(i 2fc <)) 

i=l 

= r) ■ c(k) + C • d(k) 

with ri = tb(°), C = &(?), c(k) = (tit) 2 "- 1 , and d{k) = Et=i( tr ■ c ( k ) and d ( k ) Mm thc 

recursion relations 

c(k + l) = (tit) 2 ^ 1 - 1 = (trtf- 1 (trtf =c(k)(trtf 1 

d(k + l) = J2(tr tf +1 - r = J2(tr tf + 2 " - r + (trt)° = d(k)(trtf + 1. 



i=l i=l 



Therefore 



fi(k + l) =/i(/c)(tr<) 2 ' +C- 



At this point we need a binary case distinction. Either (1.) dg(/i(/c)) is uniformly bounded by dg(C), 
implying peridocity of i, or (2.) dg(p(k)) is unbounded and passes dg(£) so no cancellation can occur 
and we can easily calculate the limit. 
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dg(p(k)) uniformly bounded by dg(£) implies dg(/i(fc)) < dg(£) for all k. Then n(t — 2 ,£) = 
dg(u(k)) is bounded and 3? = 0. But n(t,£) bounded also implies t periodic and therefore 
dg(tr t) = 0. Thus we have = = dg(tr t). 

If dg(//(/c)) is not uniformly bounded by dg(C), there exists a /z(fco) with dg(/i(fco)) > dg(C). 
Then dg(^(fc + 1)) = dg(/z(A;o)) + dg(tr t)(2 k °) and by recursion dg(^i(fc + 1)) = dg(/i(fc)) + 
dg(tr t)(2 k ), Vfc > fc . Now we can calculate the limit: 



lim l-dg(t 2k bO) = lim idg( M (fc)) 

k— >oo 2 k^-oo 2 




k — ko 



dg(/i(fc ))+ dg(trt)(2 fe - 1 ; 

, fc — fco 

& pdg(Mfco)) + fc lim ^dg(trt)^ 



=0 

= dg(trf). 

Thus in all cases = dg(tr t). □ 



Now we can proceed to the proof of Theorem IIV.8I 

Proof of Theorem \IV.8\ As shown in Theorem IIV.61 a stabilizer state of stabilizer generator length 
In + 1 encodes n maximally entangled qubits with respect to a bipartite cut. In Lemma IIV.9I we 
showed that the minimal length of a stabilizer generator grows asymptotically with 2 • dg(tr t) under 
the action of a CQCA T. Together these results prove the theorem. □ 

Figure [3] illustrates this behavior for different CQCAs. 
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FIG. 3: Entanglement generation for the stabilizer state with w(£) = 02 ® <J\ ® 02 in a bipartite setting with 
different CQCAs. One can see that entanglement can also be destroyed, but grows asymptotically linear with 
the number of timesteps n. The coefficient is given by the degree of the trace of the CSCA matrix. 



We can also calculate the entanglement of a finite region, i.e., L consecutive spins, with the rest 
of the chain. To do this calculation, we use the same method as above, and arrive at the following 
theorem. 
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Theorem IV. 10. Given a pure translation-invariant stabilizer state of stabilizer generator length 
2n + 1, a region of length L shares In maximally entangled qubit pairs with the rest of the chain if 
2n < L and L qubits pairs if 2n > L. 

Proof. The proof works exactly as in the bipartite case. In the case 2n < L the cut stabilizers are only 
cut on one side. But all stabilizers that are cut on the left side commute with those cut on the right 
side. Thus we have two independent cuts of the bipartite case and therefore 2n pairs of maximally 
entangled qubits. In the case 2n > L some stabilizers are cut on both sides. We arrange them in 
a (2n + L) x 2L-matrix like in the proof of Theorem IIV.6I and use the same technique to produce 
the mutually commuting anticommuting pairs which encode the qubits. We always find L pairs of 
maximally entangled qubits. □ 

For the evolution of entanglement under the action of a CQCA T, this means that starting with 
a product stabilizer state, the entanglement grows with 2 ■ tvt until it reaches L. Then it remains 
constant. If we start with a general translation-invariant stabilizer state, the entanglement might 
decrease at first. After some time it starts increasing and reaches L, where it remains if the CQCA is 
not periodic. Results are shown in Figure UJ 
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FIG. 4: Evolution of entanglement of a subchain of 30 consecutive spins for an initial stabilizer state with 
W (C) = o"2 ® <Ti ® ai ® <Ti ® ai ® cti <8) <T2 for different CQCA actions. The entanglement first grows as in the 
bipartite case, but then saturates at 30 qubit pairs. 



C. Entanglement generation starting from translation-invariant quasifree states 

In this section we study the entanglement generation of the glider automorphism G acting on 
a family of pure quasifree states that interpolates between the all-spins-up state (discussed in the 
previous section) and a glider-invariant state (discussed in Section UlI D 2\ . 

Let w be a pure translation-invariant quasifree state, and let wri^i denote its restriction to the 
lattice points {1,2,..., L}. The entanglement entropy of the restricted state can be calculated 

from the eigenvalues {Aj}i = i...2i, of the restricted majorana two-point matrix [-^n,?n]^ L m=1 by the 
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formula (32l.[33|: 

2L 

S , (w [li£] ) = -5^Ai/2Iog(V2)- (37) 

i=l 

The family of states that we will consider as initial states are the pure translation-invariant quasifree 
states uja described by the symbol (see Section HlID2p ): 

1 - X[-nA,0](p) + X[0,KA](P) *[1 - X[-7rA,0](p) ~ X[0,7rA](p)} 
- X[-*A,0](P) - X[0,izA](p)] 1 - X[-wA,0](p) + X[0,tt A](P) \ ' 

where X[a,b] denotes the characteristic function of the interval [a, b], and A is some real number between 
and 1. The state corresponding to A = is the all-spins- up state, while the state corresponding 
to A = 1 is a glider-invariant state. We have shown in Section IIII D 21 that by applying the glider 
automorphism n-times on lja one obtains a quasifree state belonging to the symbol 

1 - X[-7rA,o]0) + X[o,ttA](p) *[1 - X[-wA,o](p) - X[o,7tA]{p)]e 2mp 
- X[-7rA,o]0) - X[o^A]{p)Vr 2inp 1 - X[-kA,o](p) + X[o,-kA\(p) 

Using this result and Formula (|37|) . we calculated numerically the entanglement generation. The 
results are shown in Figures [5] and |U We can observe that the entanglement generation is linear in 
time, its rate is maximal when A — 0, and the rate can be arbitrarily small (when A approaches 1). 
This is illustrated in Figure [5j For longer subchains it takes more time steps for the entanglement to 
saturate. We show this in Figure [S] 



Q (ua \p) = 



200 ? 



to 
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FIG. 5: The entanglement entropy of a subchain of length 200 after applying the glider time-evolution t times. 
Different colors mark the different initial of the parameter A of the initial quasifree state. 
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FIG. 6: The entanglement entropy of a subchain of L consecutive spins in dependence of the time steps. For 
larger L more entanglement is possible. Different colors refer to different length of the chain. A is fixed at 0.9. 
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We state here some results for fractal CQCAs. All results are stated for CQCAs, because it is more 
convenient than the pure phase space formulation. Nevertheless, the results are also true for CSCAs 
and the proofs use the phase space formulation. 

Lemma A.l. A finite tensor product of only one kind of Pauli matrices (and the identity) occurs at 
most once for every Pauli matrix o~j, j = 1, 2, 3 in the history of any non-periodic CQCA T. 

Proof. We begin with the case j = 1. For two observables of the form (g) <7j, i = 1,0 to occur in the 
same time evolution of a CQCA T, the condition 



which is a periodic automaton. The case j = 3 works analogous. For j — 2 we employ the fact, that 



Appendix A: Proofs and technicalities 



1. Some results for fractal CQCA 




has to be true. It immediately follows, that 
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we can build a CQCA T which fulfills (jA.lj) via 




t = 

from any CQCA B, that fulfills &(") = (^). Moreover, as the conjugation with a CQCA does not 
change the trace, all CQCAs B have to be periodic. □ 

Lemma A. 2. Let T be a fractal CQCA on a spin chain and let §§Ui be a finite tensor product of 
Pauli matrices. For every k G N there exist an m € N such that T m (^) Ui contains at least k Pauli 
matrices. 

Proof. We assume that the number of elements is bounded by some k max . The area over which these 
kmax elements are distributed is not bounded: If it were, the elements would either be restricted to 
a finite area for an infinite number of time steps, implying periodicity, or the area would move as a 
whole implying gliders for some power of T, which is not possible as shown in Lemma III. 151 So we see 
that any group of Pauli matrices will eventually be distributed over any area. But as we require the 
number of elements to be bounded, the distance between any two groups of Pauli matrices becomes 
larger than the neighborhood of the automaton. Then the starting argument applies to each of the 
new groups and forces them to break apart further until only isolated single-cell observables are left. 
But these will expand, thus the number of elements can't be bounded. □ 



2. Proof of theorem IIV.6I 



Proof of Theorem \IV.6\ We use the criterion of Theorem lIV.4l and explicitly construct the pairs w(£j), 
w(t^) using methods from stabilizer codes for quantum error correction [301 ] . As said in Section HVBI 
only stabilizer generators localized on both sides of the cut are elements of the correlation group Sab- 
The projections of all other stabilizer generators are just the stabilizer generators themselves which 
trivially commute with all other stabilizer generators and their projections on A resp. B. We now use 
the phase space representation of the Pauli products and build the following 2n x 4n-matrix from the 
cut translates of the stabilizer generators: 



(c 



t-n+l 



V 



o 



e + - 



o c n 
o 



£71 — 2 



o c n J 



Let us assume that the outermost element is a <7i 14 . Then £7" = 1 and £l™ = 0. From Section TlVBI 
we know, that at least one £!i 7^ 0. Let the i-th diagonal of the right part be the first non-zero one. 
We also use the reflection invariance of £ to replace by £ 3 and get the following matrix: 



/ 1 ^r 1 





Vo 



1 c+ 



tn-2 



1 a~ 
1 











\ 



C 



/ 



The other cases work equivalently. 
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Now we can perform the Gaussian algorithm on the matrix to obtain an identity matrix in the left 
part. As the rows are shifted copies of the first row, all operations will also be applied in a shifted 
copy. If we would add the third row to the first, we would also add the fourth to the second and so 
forth. We only add rows to rows above, because the lower left part of the matrix is already zero. We 
therefore get 



( I ••• 
1 ••• 







1 
1 



o ••• o c 



o 



... £7.-1 ^ 
n+i 



o J 



The i-th row of the right part of the matrix remains unchanged, so ( n+t 
get operators of the form 



1. Wc therefore 



-n+i— 1 



Sk = w(£ fe ) = a\ 



(g) V 



-n-\-i 



1< k < 2n, 



\j = -n+l 



\j = -n+i+l 



where the cr? can only be 03 or gq = 1. Wc can easily see, that {s^, Sk+i-i} = 0. As i < n 
we can always find n pairs of anticommuting operators. Unfortunately these pairs do not necessarily 
commute with other pairs. But through multiplication of operators we can find new pairs, which fulfill 
the necessary commutation relations. To show this we create a (symmetric) matrix Ci j — cr(si, Sj) of 
commutation relations. 



C 





Sl 


Si 


s 3 


Sj+i-1 ■ ■ ■ 


Sl 





1 


? 


? 


St 


1 





? 


? 


S J 


? 


? 





1 


Sj+i-1 


? 


? 


1 






A "1" stands for anticommutation, a "0" for commutation. In the end we want all operators from 
different pairs to commute, so all positions denoted by question marks should get a zero entry. If we 
multiply operators, the corresponding rows and columns are added. Through these operations we can 
bring the commutation matrix to the form 



C 





Sl 


Si 


Sj 


• • • 


Sl 





1 








Si 


1 











Sj 











1 


Sj+i-1 








1 






To show that this is possible, we will consider a prototype of such an operation. Given the matrix 



C 





Sl 


Si 


S 3 


Sj+i-1 • • • 







1 


Cl3 


C14 


Si 


1 





C23 


C 2 4 


Sj 


Cl3 


C23 





1 


Sj+i-1 


C14 


C24 


1 
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We now pick a nonzero and do the following. If Cjj = 1 and i odd, we add row i + 1 to row j 
and the same for the columns. If c^- = 1 and i even, we add row i — 1 to row j and the same for the 
columns. This only changes the one dj we are considering, the others remain unchanged. After each 
step we get a new matrix C and pick another nonzero Cjj from the same 2x2 block (in this example 
we only have one block). By doing this for all blocks in the first two rows, we create pairs of operators 
that commute with the first pair. Now we have to check if the process destroyed the anticommutation 
within the pairs. The diagonal entries of the matrix trivially stay zero, because all operators commute 
with themselves. We only have to check the other elements of the block (due to the symmetry, we 
only have to check one). So if C13 = 1 we get 1^1 + c 2 4. We can write 1 1 — >■ 1 + C13C24. Including the 
whole block of dj we get 1^1 + ci 3 c 2 4 + C14C23 + C23C14 + C24C13 = 1 as all operations are carried out 
modulo 2. The new pairs thus fulfill the anticommutation condition. We can repeat this process for 
the new pairs until all operators from different pairs commute. We started with 2n operators, thus 
we arrived at n pairs which together with their counterparts on the other subsystem encode n pairs 
of maximally entangled qubits. □ 
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